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Abstract— Based on published data, we have collected information about Galactic 
\q \ maser sources with measured distances. In particular, 44 Galactic maser sources located in 

star-forming regions have trigonometric parallaxes, proper motions, and radial velocities. 
In addition, ten more radio sources with incomplete information are known, but their 
parallaxes have been measured with a high accuracy. For all 54 sources, we have calculated 
the corrections for the well-known Lutz-Kelker bias. Based on a sample of 44 sources, 
we have refined the parameters of the Galactic rotation curve. Thus, at R = 8 kpc, 
the peculiar velocity components for the Sun are (Uq,V q ,Wq) = (7.5,17.6,8.4) ± 
(1.2, 1.2, 1.2) km s -1 and the angular velocity components are uj = —28.7 ± 0.5 km 
s" 1 kpc" 1 , u' Q = +4.17 ± 0.10 km s -1 kpc~ 2 , and u$ = -0.87 ± 0.06 km s -1 kpc" 3 . The 
corresponding Oort constants are A = 16.7 ± 0.6 km s _1 kpc -1 and B = —12.0 ± 1.0 km 
s _1 kpc -1 ; the circular rotation velocity of the solar neighborhood around the Galactic 
center is Vn = 230 ± 16 km s _1 . We have found that the corrections for the Lutz-Kelker 
ly-j ■ bias affect the determination of the angular velocity loq most strongly; their effect on the 

remaining parameters is statistically insignificant. Within the model of a two-armed spiral 
pattern, we have determined the pattern pitch angle i = — 6.°5 and the phase of the Sun 
^ \ in the spiral wave Xo — 150°. 



1 INTRODUCTION 



X ■ 

The observed velocities of various young Galactic objects (HI clouds, OB stars, young open 
star clusters) are used to study the kinematics of the Galaxy and the spiral structure of 
the Galactic disk (Bobylev et al. 2008). 

For this task, Galactic maser sources with trigonometric parallaxes measured by VLBI 
are extremely important (Reid et al. 2009a; McMillan and Binney 2010; Bobylev and 
Bajkova 2010; Stepanishchev and Bobylev 2010). We have in mind the maser sources 
associated with the youngest Galactic stellar objects (either protostellar objects of various 
masses, or very massive supergiants, or T Tauri stars). The number of measured parallaxes 
for such masers already exceeds 50 and with a complete data set (with proper motions 
and radial velocities) is 44. 

Studying their distribution and kinematics with allowance made for the distance errors 
is of great importance. First of all, it is necessary to take into account the Lutz-Kelker 
(1973) bias. The corrections for this bias usually play an important role in determining 
the trigonometric parallaxes by optical methods, when the reference objects are distant 
stars rather than extragalactic sources. In this procedure for absolutizing the 
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measured relative parallaxes is needed. For this purpose, the corrections calculated from 
the Lutz-Kelker bias are usually applied. 

In the case of Galactic maser sources, the situation is different. Their trigonometric 
parallaxes are determined by VLBI with referencing directly to distant quasars (two or 
three quasars are commonly used). Therefore, such parallaxes are immediately obtained 
as absolute ones. 

At the same time, the VLBI parallaxes of masers were measured with a certain random 
error that is, on average, very small, about 5%. However, there are exceptions and, in 
several cases, the error exceeds 20%. As is well known (Lutz and Kelker 1973), the 
knowledge of the true parallax depends strongly on the relative error of the measured 
parallax. It can be surmised that the parameters of the Galactic rotation curve and the 
parameters of the spiral density wave determined from maser sources can be distorted 
because of the uncertainties in their distances. 

Great experience in taking into account the Lutz-Kelker bias has been gained. We 
will point out only some of the papers. Hanson (1979) applied the method in practice to 
calibrate the luminosity based on trigonometric parallaxes. In the literature, the method 
is occasionally called the Lutz-Kelker-Hanson method. Maiz-Apellaniz (2001) determined 
such corrections for 0-B5 stars with trigonometric parallaxes from Hipparcos (1997), 
with the Sun's height above the Galactic plane having been taken into account. Verbiest 
et al. (2010) proposed a modification of the method that took into account the luminosity 
function as applied to 57 Galactic pulsars with measured trigonometric parallaxes. Note 
that the Lutz-Kelker bias is a statistical effect. A very accurate model of the expected 
distribution of stars in the Galaxy is required to properly take the Lutz-Kelker bias into 
account. 

The goal of this paper is to calculate the corrections for the Lutz-Kelker bias for masers 
with measured trigonometric parallaxes and to study the effect of these corrections on the 
parameters of the Galactic rotation curve and the parameters of the spiral density wave 
being determined. 



2 METHOD 

2.1 Uniform Distribution in Space 

In its classical form, the Lutz-Kelker bias refers to a uniform distribution of stars in an 
infinite space (Lutz and Kelker 1973). We assume that the errors in the distribution of 
the measured parallax relative to its true value are distributed normally: 

Here, 7r = 3.14..., a y is the measurement error of the parallax, wq and w are the 
measured and true parallaxes, respectively. We are interested in the distribution of true 
parallaxes w if a sufficient number of stars with a certain measured value of wq are 
available. If the stars are uniformly distributed in space, then their number in an interval 

(r, r + dr) is 

N(r)dr = 4vrr 2 dr, (2) 
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or 

, , , 4irdm , . 

N(w)dw = — . (3) 



According to the Bayes theorem, 



P{w\Wq) = ? r , 4 

where p(cc7o) may be set equal to a constant. We find that 

pM^ ) = -exp^ 2^ J " ( 5 ) 

The constant k is chosen in such a way that the normalization condition is fulfilled: 

p{w\wQ)dw = 1. (6) 

The lower limit vj m i n implies the maximum distance at which a given star can be observed 
in principle from the considerations of the stellar system's size and the telescope's limiting 
magnitude. In the case of integration from zero, the integral diverges. This reflects an 
infinite number of stars in an infinite space for the chosen uniform distribution. 

The parallax corrected for the Lutz-Kelker bias is the expectation of distribution (151) : 

oo 

ttcorr = / VJp{w\Wo)(\vj . (7) 



2.2 An Exponential Disk with an Observer at the Periphery 

Our objective is to study the objects distributed in the Galactic disk. We assume an 
exponential density distribution along the Galactic radius, with the distribution along 
the vertical Z axis following the law of the hyperbolic secant squared: 

p(#)=p exp(-|)ch- 2 (J^, /i = 3kpc, Z = 0.25 kpc, (8) 

where R is the Galactocentric distance, h and Zq are the radial and vertical scale lengths 
of the disk, respectively; the constant po, the density at the center, is a scale factor 
and its value is of no importance for our problem. The Sun is at Rq = 8 kpc from the 
Galactic center. When observing from the Sun along the Galactic longitude /, the density 
distribution takes the form 



(rcos/ - i? ) 2 + r 2 sin A 
p(r, I) = po exp -- * I. (9) 

The parallax corrected for the Lutz-Kelker bias is calculated from (J7J) with the distribution 
of true parallaxes at given observables: 

p( WWo ) oc ^^exp . (10) 

ro 4 \ 2 o I 
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Ta6jiHU,a 1: Data on the sources that were not included in the "kinematic" sample 



Source 


ot 


8 


^ y~"UJ J 


n* (u ,, ) 


LLx((7,, ) 


V„ ( (Tt/ \ 


Ref 


G 23.66-0.13 


278.7149 


-8.3059 


.313(.039) 


-1.32(.02) 


-2.96(.03) 


83(3) 


(a) 


G 9.62+0.20 


271.5611 


-20.5255 


.194(.023) 


-.58(.05) 


-2.49(.27) 


0.1(2) 


(b) 


Sgr B2N 


266.8330 


-28.3720 


.128(.015) 


-.32(.05) 


-4.69(.ll) 


64(5) 


(c) 


Sgr B2M 


266.8340 


-28.3845 


.130(.012) 


-1.23(.04) 


-3.84(.ll) 


61(5) 


(c) 


IRAS 20126+41 


303.6084 


41.2257 


.61(.02) 


-2.0(.l) 


l.O(.l) 


-3.5(4) 


(d) 


G 48.61+0.02 


290.1299 


13.9237 


.199(.007) 


-2.76(.04) 


-5.28(.ll) 


19(1) 


(e) 


MSXDC G034.4+0 


283.3292 


1.4022 


.643(.049) 


-.25(.18) 


.00(.18) 


57(5) 


(0 


EC 95 


277.4912 


1.2128 


2.41(.02) 


.70(.02) 


-3.64(.10) 




(g) 


IRAS 05137+3919 


79.3073 


39.3722 


.086(.027) 


.30(.10) 


-.89(.27) 


-26(3) 


(h) 


G 27.36-0.16 


280.4627 


-5.0287 


.125(.042) 


-1.81(.08) 


-4.11(.26) 


92.2(5) 


(i) 



Note, a and 5 are given in degrees, the parallax w is in mas, /x* = fi a cos 8 and are in mas 
yr _1 , the radial velocity V r = V r {LSR) is in km s -1 ; the letters mark the references to papers: 
(a) Bartkiewicz et al. (2008); (b) Sanna et al. (2009); (c) Reid et al. (2009); (d) Moscadelli 
et al. (2011); (e) Nagayama et al. (2011b); (f) Kurayama et al. (2011); (g) Dzib et al. (2010); 
(h) Honma et al. (2011); (i) Xu et al. (2011). 



3 DATA 

The trigonometric parallaxes and proper motions of Galactic maser sources are determined 
by several research groups using long-term VLBI observations. These include the Japanese 
VERA (VLBI Exploration of Radio Astrometry) project on the observation of Galactic 
H20 maser sources at 22 GHz and SiO masers (these are very few among the young 
objects) at 43 GHz. Note that the higher the frequency, the higher the resolution and 
the more accurate the observations. Methanol (CH30H) masers are observed at 12 GHz 
using the American NRAO VLB A and at 6.7 GHz using the European VLBI Network. 
At present, Galactic masers are observed by these groups of researchers within the unified 
BeSSeL project (Brunthaler et al. 2011) aimed at studying the Galactic structure. VLBI 
observations of a number of young radio stars in continuum at 8.4 GHz (Dzib et al. 2010) 
are also conducted to study the structure and kinematics of the solar neighborhood (the 
Gould Belt, molecular clouds). 

Bobylev and Bajkova (2010) and Stepanishchev and Bobylev (2011) analyzed a sample 
of 28 masers with measured trigonometric parallaxes drawn from published data. By 
now, the amount of such data has increased considerably— about 20 new measurements 
in various star-forming regions have been published. First of all, we use the input data on 
44 masers with measured trigonometric parallaxes, proper motions, and radial velocities 
from Bajkova and Bobylev (2012), where corresponding references and explanations are 
given. We arbitrarily call this sample "kinematic", because on its basis we determine the 
Galactic rotation parameters. 

Table 1 gives the parameters of the sources most of which have small relative errors 
in the parallaxes, but, for various reasons, they were not included in the "kinematic" 
sample. Two sources are almost at the Galactic center, Sgr B2N and Sgr B2M (Read 
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Phc. 1: Projections of masers onto the Galactic plane (XY) from the observed data 
(crosses) and corrected for the Lutz-Kelker bias (circles); the large filled square marks the 
position of the Sun; the Galactic Center is at the coordinate origin. 



et al. 2009a). Other sources have large (>40 km s r ) residual velocities: these include 
G23.66.0.13 (Bartkiewicz et al. 2008), MSXDC G034.43+00.24 (Kurayama et al. 2011), 
G48.61+0.02 (Nagayama et al. 2011), and IRAS 20126+4104 (Moscadelli et al. 2011). 
Since the star-forming region G9. 62+0. 20 (Sanna et al. 2009) is closely connected with 
the 3-kpc spiral arm, it has a large peculiar Galactocentric radial velocity, |Vr| ~ 50 
km s _1 . In addition, we included the radio star EC 95 (Dzib et al. 2010) in our sample, 
for which there are no radial velocity measurements as yet. Its VLBI observations were 
carried out in continuum at 8.42 GHz. It is hypothesized that this is a binary system 
consisting of a Herbig-Haro protostar with a mass of about 4 — 5 M Q and a T Tauri 
low-mass companion. 

The stars from Table 1, especially those farthest from the Sun, are of interest in 
studying the effect of the Lutz-Kelker bias on the determination of the spiral-structure 
parameters. 

4 RESULTS AND DISCUSSION 



The observed (w &s) an d Lutz-Kelker-biascorrected (w corr ) parallaxes for all 54 maser 
sources are given in Table 2; their distribution in projection onto the Galactic XY plane 
is shown in Figs. 1 and 2. As can be seen from Table 2, for most of the maser sources, the 
calculated Lutz-Kelker bias corrections are negative (the heliocentric distances to them 
should be increased). The correction is largest for IRAS 16293—2422, but this source is 
one of the closest to the Sun. Therefore, the correction found has no critical effect on 
the determination of the (distance-dependent) parameters under consideration. As can 
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Phc. 2: Two-armed spiral pattern constructed from masers in the Perseus and Outer Arms 
(thick circles); the remaining objects from our list are indicated by the crosses. 



be seen from Fig. 1, a significant change in the distance is achieved for several sources 
near the Galactic center. The source farthest from the Sun, IRAS 05137+39, is greatest 
interest. The relative error in the parallax for it is ?^30%; since the source is close to the 
direction of the Galactic anticenter, the correction for the Lutz-Kelker bias is positive (its 
heliocentric distance should be reduced). 

Figure 1 clearly illustrates the essence of the method. More specifically, the existence 
of small measurement errors in the parallaxes gives grounds to expect that the most 
probable distance to a star should be closer to the region with the highest density of stars 
(in our case, this is the Galactic center region). 

4.1 Spiral Pattern 

The parameters of the spiral pattern were determined by the method described in 
Stepanishchev and Bobylev (2011). The essence of the method consists in the following. 
We proceed from the well-known equation (Lin and Shu 1964) 

X = m(ctgi-hi{R/R }-6)+Xo, (11) 

where x is the phase of the object in the spiral density wave, xo is the phase of the Sun 
in the spiral wave, R is the distance from the object to the Galactic rotation axis, 9 is the 
Galactocentric longitude of the object, and m is the number of spiral arms (here, we take 
m = 2). Having chosen two or three groups of objects located along the arms, we estimate 
the most probable phase of the Sun and the pitch angle of the spiral pattern defined by 
Eq. (11). For this purpose, we search for the minimum of the functional 

C 2 (*,xo) = £(to) 2 > (12) 
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where Xj is the phase of the jth object reduced to the range (— 7r/2,7r/2) and pj is the 
weight that is chosen based on the number of objects in a given arm in such a way that all 
arms have the same weight. Here, the masers are assumed to be theoretically concentrated 
to the spiral arm centers, i.e., the sum of the squares of the objects' phases in the density 
wave is minimized by taking into account the fact that the phase x = 0°. corresponds to 
the arm center. 

Owing to the increase in the amount of observational data, these parameters can be 
calculated with a higher accuracy. Thus, three masers are already known in the Outer 
Arm against two in the previous paper, while eleven masers are known in the Perseus 
Arm against seven. In Fig. 2, the objects involved in choosing the model parameters are 
indicated by the filled circles and the remaining ones are indicated by the crosses. The 
pitch angle and the Sun's phase were found to be, respectively, i = —6°. 5 and xo = 150°. 

Since the Lutz-Kelker bias corrections for most of the masers belonging to the spiral 
arms (from which the parameters of the spiral pattern were determined) are insignificant, 
they were disregarded in the case under consideration. Note once again the source IRAS 
05137+39, which is behind the Outer Arm (at ^20 kpc from the Galactic center). It was 
not used in choosing the model parameters; nevertheless, it falls nicely on the logarithmic 
spiral fitted into the Perseus and Outer Arms (Fig. 1 and Fig. 2). In their kinematic 
analysis of a sample of Galactic masers, Bajkova and Bobylev (2012) found a slightly 
smaller pitch angle of the spiral wave, i = — 5lo!g° (f° r m = 2), and a phase of the Sun 
in the spiral wave fairly close to our result, x® — — 147lf 7 o (measured from the Carina- 
Sagittarius Arm). At the same time, the spiral pitch angle we derived differs significantly 
from the result of analyzing the samples of masers in the Perseus and Outer Arms obtained 
by Sakai et al. (2012) based on a different method (with a separate analysis of each arm). 
These authors found % = —17°. 8 ± 1°.7 for the Perseus Arm and % = —11°. 6 for the Outer 
Arm. This discrepancy in determining the pitch angle can be explained by the fact that 
the results of Sakai et al. (2012) more closely correspond to the four-armed model of the 
spiral pattern (in our case, we will then have i = —13°). On the whole, we can conclude 
that considerably larger statistics is required for a reliable determination of the spiral 
pattern parameters. 

4.2 Galactic Rotation Parameters 

First, we redetermined the parameters of the Galactic rotation curve using data on 44 
masers. Our technique of its determination based on a Taylor expansion of the angular 
velocity in terms of the Galactocentric distance was described in detail in Bobylev and 
Bajkova (2010) and Bobylev and Stepanishchev (2011). At a fixed Galactocentric distance 
of the Sun, R = 8 kpc (with an uncertainty of 0.5 kpc), the peculiar velocity components 
of the Sun (U Q ,V Q ,W e ) = (7.5,17.6,8.4) ± (1.2,1.2,1.2) km s^ 1 and the following 
parameters of the Galactic rotation curve were obtained: 

oj = -28.7 ±0.5 km s" 1 kpc^ 1 , 

u' = +4.17 ±0.10 km s' 1 kpc" 2 , (13) 
u" = -0.87 ± 0.06 km s" 1 kpc" 3 , 
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where cu is the angular velocity of Galactic rotation at R = R , u' and u>q are its 
corresponding derivatives. The Oort constants are A = 16.7 ± 0.6 km s _1 kpc -1 and 
B = —12.0 ± 1.0 km s -1 kpc -1 ; the circular rotation velocity of the solar neighborhood 
around the Galactic center is Vq = 230 ± 16 km s _1 . 

Note that Stepanishchev and Bobylev (2011) found the analogous parameters using 
only 28 masers to be the following: (U e , V e , W Q ) = (8.5, 17.1, 8.9) ± (1.6, 1.6, 1.6) km s~\ 
and cuo = -30.4 ± 0.7 km s" 1 kpc" 1 , u' Q = 4.23 ± 0.13 km s" 1 kpc" 2 , u% = -1.01 ± 
0.06 km s -1 kpc~ 3 . It can be seen that solution (13) is more reliable, because, as would be 
expected, the random errors of all the parameters being determined decreased. Solution 
(13) virtually coincides with results of analyzing this sample of masers from Bajkova and 
Bobylev (2012), where a detailed comparison of these values with the results of other 
authors can be found. 

After applying the Lutz-Kelker correction, we obtained the following parameters of 
the rotation curve: (U Q , V Q , W & ) = (6.6, 17.9, 8.5) ± (1.4, 1.4, 1.4) km s" 1 and 

u = -27.5 ±0.5 km s -1 kpc" 1 , 

uj' = +3.98 ±0.11 km s" 1 kpc" 2 , (14) 
u" = -0.83 ± 0.06 km s" 1 kpc" 3 , 

Hence we see that after applying the correction, the angular velocity w changed by about 
2.5(7, while its derivatives and the Sun's peculiar velocity components change within the 
error limits (i.e., insignificantly). 

5 CONCLUSIONS 

Based on published data, we collected information about Galactic maser sources with 
well measured distances. In particular, 44 Galactic maser sources located in star-forming 
regions have a complete set of kinematic observational data-the VLBI trigonometric 
parallaxes and proper motions; the mean radial velocities derived from CO observations 
are also known for them. In addition, ten more radio sources with incomplete information 
are known, but their parallaxes were measured with a sufficiently high accuracy. 

We calculated the corrections for the well-known Lutz-Kelker bias for all 54 sources. 
For this purpose, we used a model with an exponential star density distribution in the 
Galactic disk. For several sources farthest from the Sun, these corrections were shown to 
be appreciable. Using the calculated corrections can be useful, for example, in studying 
the membership of masers in the Galactic bar or specific spiral arms, etc. 

Based on a sample of 44 sources, we refined the parameters of the Galactic rotation 
curve. At R0 = 8kpc, we determined the peculiar velocity components of the Sun 
(Uq,Vq,W q ) = (7.5,17.6,8.4) ± (1.2,1.2,1.2) km s" 1 and the following parameters: 
coo = -28.7 ±0.5 km s" 1 kpc -1 , cu' = ±4.17 ±0.10 km s" 1 kpc~ 2 h u% = -0.87 ±0.06 km 
s _1 kpc~ 3 . The corresponding Oort constants are A = 16.7 ± 0.6 km s _1 kpc" 1 and 
B = —12.0 ± 1.0 km s _1 kpc -1 ; the circular rotation velocity of the solar neighborhood 
around the Galactic center is V = 230 ±16 km s -1 . We found that when using our sample, 
the corrections for the Lutz-Kelker bias affect the determination of the angular velocity 
of Galactic rotation uq (at the 2.5<r level), while these corrections have no statistically 
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significant effect on the other model parameters being determined. This is primarily due 
to the high accuracy of the measured trigonometric parallaxes for the masers used. Within 
the model of a two- armed spiral pattern, we determined the pattern pitch angle i = —6°. 5 
and the phase of the Sun in the spiral wave xo — 150°. 
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Ta6jiiiiia 2: Observed and Lutz-Kelker-bias-corrected parallaxes (in 
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Table 2 (the end) 



Source 


^obs 


^ 1 COTT 


Correction 


HW9 CepA 


1.430 (0.070) 


1.417 


-0.013 


Sgr B2N 


0.128 (0.015) 


0.122 


-0.006 


Sgr B2M 


0.130 (0.012) 


0.126 


-0.004 


IRAS 05137+3919 


0.086 (0.027) 


0.093 


0.007 


G 27.36-0.16 


0.125 (0.042) 


0.103 


-0.022 


G 48.61+0.02 


0.199 (0.007) 


0.198 


-0.001 


IRAS 20126+4104 


0.610 (0.020) 


0.608 


-0.002 


G 9.62+0.19 


0.194 (0.023) 


0.177 


-0.017 


MSXDC G034.43+00.24 


0.643 (0.049) 


0.626 


-0.017 


G 23.66-0.13 


0.313 (0.039) 


0.286 


-0.027 


EC 95 


2.410 (0.020) 


2.409 


-0.001 
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